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ABSTRACT 

The unknown nature of “dark energy” motivates continued cosmological tests of large- 
scale gravitational physics. We present a new consistency check based on the relative 
amplitude of non-relativistic galaxy peculiar motions, measured via redshift-space 
distortion, and the relativistic deflection of light by those same galaxies traced by 
galaxy-galaxy lensing. We take advantage of the latest generation of deep, overlapping 
imaging and spectroscopic datasets, combining the Red Cluster Sequence Lensing Sur¬ 
vey (RCSLenS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), 
the WiggleZ Dark Energy Survey and the Baryon Oscillation Spectroscopic Survey 
(BOSS). We quantify the results using the “gravitational slip” statistic Eq, which we 
estimate as 0.48 ±0.10 at 0 = 0.32 and 0.30 ±0.07 at z = 0.57, the latter constituting 
the highest redshift at which this quantity has been determined. These measurements 
are consistent with the predictions of General Relativity, for a perturbed Friedmann- 
Robertson-Walker metric in a Universe dominated by a cosmological constant, which 
are Eg = 0.41 and 0.36 at these respective redshifts. The combination of redshift-space 
distortion and gravitational lensing data from current and future galaxy surveys will 
offer increasingly stringent tests of fundamental cosmology. 
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1 INTRODUCTION 

A wide set of cosmological observations suggest that the dy¬ 
namics of the Universe are currently dominated by some 
form of “dark energy”, which in standard Friedmann- 
Robertson-Walker (FRW) models is propelling an acceler¬ 
ation in late-time cosmic expansion (e.g. Planck collabora¬ 
tion 2015a, Aubourg et al. 2015, Betoule et al. 2014). The 
physical nature of dark energy is not yet understood, and a 
widely-considered possibility is that the nature of gravita¬ 
tion differs on large cosmological scales from the predictions 
of General Relativity (GR) in an FRW metric. As a result, 
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a key task for current cosmological surveys is to construct 
observations to test for such departures. 

Gravitational physics produces a rich variety of observ¬ 
able signatures that can be used for this purpose. The most 
precisely-measured signal results from the “peculiar mo¬ 
tions” of galaxies as they fall toward overdense regions as 
non-relativistic test particles in a perturbed FRW metric. 
These motions produce correlated Doppler shifts in galaxy 
redshifts that manifest themselves as an overall anisotropy 
in the measured clustering signal as a function of the angle to 
the line-of-sight (Kaiser 1987), known as redshift-space dis¬ 
tortion (RSD). The amplitude of this anisotropy has been 
accurately measured in a number of galaxy surveys across 
a range of redshifts and allows the growth rate / of cosmic 
structure, which describes the gravitational amplification of 
density perturbations, to be inferred. To date, these mea- 


© 0000 RAS 


2 Blake et al. 


surements are broadly consistent with the prediction of the 
standard cosmological model (e.g. Blake et al. 2011, Beutler 
et al. 2012, Reid et al. 2012, de la Torre et al. 2013, Samushia 
et al. 2014). 

A highly complementary route for probing gravitational 
effects is to study the deflections of relativistic test par¬ 
ticles such as photons, which are additionally sensitive to 
the curvature of space produced by density perturbations. 
These deflections, known as gravitational lensing, may be 
measured in a statistical sense using the correlations im¬ 
printed in the apparent shapes of background galaxies be¬ 
hind foreground lenses in deep imaging surveys. The level 
of the signal is determined by both the amplitude of the 
density fluctuations around the lenses (again reflecting the 
growth of structure with redshift) and the relative distances 
of the source-lens systems, both of which may be predicted 
by a given cosmological model (for a review, see Bartelmann 
& Schneider 2001). Recent projects such as the Canada- 
France-Hawaii Telescope Lensing Survey (CFHTLenS) have 
allowed a suite of such tests to be carried out (e.g. Heymans 
et al. 2013, Simpson et al. 2013) by providing deep, wide, 
high-resolution imaging. 

Zhang et al. (2007) proposed that a powerful gravita¬ 
tional consistency check might be performed by using the 
same set of galaxies to trace non-relativistic gravitationally- 
driven motion using redshift-space distortion, and to serve 
as foreground lenses for probing the relativistic deflection 
of light from background sources. In this way, it could be 
ascertained whether the relative amplitude of these two ef¬ 
fects, driven by the same underlying density perturbations 
traced by the lenses, was consistent with the prediction of 
GR assuming a perturbed FRW metric for a given set of 
cosmological parameters including the matter density fl m . 
This can be achieved by measuring a quantity known as the 
“gravitational slip” Eg(R) as a function of physical scale R, 
which is constructed from the galaxy-galaxy lensing signal 
and the redshift-space distortion and clustering amplitude 
of the lenses. Standard perturbed GR cosmology predicts 
that a scale-independent value Eg = Qm/f should be re¬ 
covered. Failure of this cross-check would indicate either the 
breakdown of linear perturbation theory, an inconsistency 
in the assumed cosmological parameters such as the matter 
density or curvature, or that a large-scale modification in 
gravitational physics was required. 

A requirement for carrying out this programme is the 
availability of galaxy spectroscopic redshift surveys and deep 
gravitational lensing imaging surveys which overlap on the 
sky. Reyes et al. (2010) performed this test using data from 
the Sloan Digital Sky Survey (SDSS) covering « 5000 deg 2 , 
using a sample of a: 70,000 Luminous Red Galaxy lenses 
at z = 0.32 and shape measurements with a surface den¬ 
sity « 1 arcmhY 2 and produced a gravitational slip mea¬ 
surement consistent with GR in a A-dominated Universe. 
The goal of our current study is to use new, significantly 
deeper, spectroscopic and imaging survey datasets to ex¬ 
tend these tests to higher source densities and more distant 
redshifts z ~ 0.6. In particular we use imaging data from 
CFHTLenS (Heymans et al. 2012) and the 2nd Red Cluster 
Sequence Lensing Survey (Gilbank et al. 2011, Hildebrandt 
et al. 2015), and overlapping spectroscopic data from the 
WiggleZ Dark Energy Survey (Drinkwater et al. 2010) and 
Baryon Oscillation Spectroscopic Survey (BOSS, Eisenstein 


et al. 2011) to carry out this investigation. We note that a 
measurement of Eg using lensing of the Cosmic Microwave 
Background was recently presented by Pullen et al. (2015). 

The structure of our paper is as follows: in Section [5] 
we introduce the theory underpinning galaxy-galaxy lens¬ 
ing and RSD in a perturbed FRW metric, the methodology 
for suppressing small-scale information which is difficult to 
model, and our test statistic for gravitational physics, Eg- 
In Section [3] we summarize our input datasets, and in Sec¬ 
tion 5 ] we present the galaxy-galaxy lensing measurements 
expressed as both the average tangential shear as a function 
of angular separation 7 t(9) and the differential surface den¬ 
sity as a function of projected physical separation AE (R), 
together with a series of tests for shape-measurement sys- 
tematics. We pay particular attention to how the full redshift 
probability distribution of each source, determined from 
photometric-redshift estimation, can be included in an unbi¬ 
ased estimator of AE (R) from source-lens pairs. In Section 
[5] we describe a large set of new N-body simulations, in¬ 
cluding self-consistent gravitational lensing, that we use for 
determining errors in our measurements and testing models. 
Finally in Section [U] we present the first cosmological impli¬ 
cations, including new measurements of Eg up to redshift 
z « 0.6. We conclude in Section [7] 

This paper is the first in a series which will dissect the 
cosmological information in overlapping deep galaxy lens¬ 
ing and spectroscopic datasets. Future studies will present 
full cosmological fits including Cosmic Microwave Back¬ 
ground data in a range of modified gravity scenarios, and the 
use of photometric-redshift and spectroscopic-redshift cross¬ 
correlations to determine simultaneously the source redshift 
distributions and cosmological parameters. 


2 THEORY: GALAXY-GALAXY LENSING 
AND REDSHIFT-SPACE DISTORTION 

2.1 Galaxy-galaxy lensing: tangential shear 

Galaxy-galaxy lensing describes the shear-density correla¬ 
tion imprinted between the shapes of background source 
galaxies and the foreground large-scale structure traced 
by the lens galaxies; it is measured by computing the 
azimuthally-averaged tangential shear of the sources as a 
function of their angular distance 9 from the lenses: 

<7 tie)) = {S g (x ) 7 t(x + 6)), ( 1 ) 

where 7 1 denotes the tangential shear component of the 
source with respect to the separation vector 9 connecting it 
to the lens, and S g (x) is the overdensity of the lens galaxies 
at position x. Fourier transforming the variables in Equation 
|T]we find that (e.g. Hu & Jain 2004) 

{lt{9)) = hJ 0 dUC ^)Mm, ( 2 ) 

where J 2 denotes the 2nd-order Bessel function of the first 
kind, and the galaxy-convergence cross-power spectrum C gK 
is given by (e.g. Bartelmann & Schneider 2001, Guzik & 
Seljak 2001, Joachimi & Bridle 2010) 

^ 312 m R 2 H ,.(! + *)„ (I ^..w^ 

C gK (£) — ^ 2 J cL, J pi(z)W(z), (3) 
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where 


W(z) 


dz' Ps(z') 


X(z') ~ X(z) 

X(z') 


(4) 


In these equations, x( z ) is the co-moving radial coordinate 
at redshift z assuming a spatially-flat Universe, P g m(k,z) 
is the galaxy-mass cross-power spectrum at wavenumber k, 
\ps(z),pi(z)] are the redshift probability distributions of the 
source and lens samples, normalized such that p(z ) dz = 
1 , and Ho and c are Hubble’s constant and the speed of 
light, respectively. These relations make the approximations 
of using the Limber equation and neglecting additional ef¬ 
fects such as cosmic magnification and intrinsic alignments. 

These equations make clear that the tangential shear is 
sensitive to the cosmological model through the distance- 
redshift relation describing the relative geometry of the 
source-lens-observer systems, and through the clustering of 
the matter overdensities traced by the lenses. Our method 
of estimating (7 t(9)) from the data is described in Section 

IXT1 


2.2 Galaxy-galaxy lensing: differential surface 
density 


In this study we are interested in “de-projecting” the angu¬ 
lar tangential shear profile in order to recover the statistics 
of the projected mass surface density E around the lenses 
as a function of projected spatial separation R. Assuming 
an axisymmetric lens, and a single lens-source pair with re¬ 
spective redshifts ( zi,z s ), the average tangential shear at 
projected separation R from the lens is given by 


< 7 * m 


A X(R, Zl ) 
E c (zz, z s ) 


(5) 


(e.g. Bartelmann & Schneider 2001), where 6 = R/\(zi) and 
the differential projected surface mass density AS is defined 
in terms of E(i?) by 


AS (R) = E(< R) - E (R), 


(6) 


where E(< R) is the average mass density within a circle of 
radius R , 


_ o f R 

E( <R) = —J R'X(R')dR'. (7) 

In Equation [5] the co-moving critical surface mass density 
E c is given by 


S ^') = 4^ 


X(z 


[x{z')-x(z)\ x(z){l + z) 


(8) 


where G is the gravitational constant and z' > z. The differ¬ 
ential projected surface density around the lenses is related 
to the lens galaxy-matter cross-correlation function £ gm (r) 
by 


AE (R) = p^ / tgmiy/lP+IP) dH , 


(9) 


where pTf is the mean cosmological matter density, related to 
the critical density p c = 3Ho/8nG = 2.77518 x 10 11 /i 2 Mq 
M pc -3 [with h = Ho/ (100 km s -1 Mpc -1 )] by pTT = pc Xl m , 
and n denotes the co-moving line-of-sight separation. 


By substituting Equations [3] and [4] into Equation [2] and 
comparing with the result of substituting Equations [8] and 
[9] into Equation [5] we can determine after some algebra how 
Equation [5] generalizes for broad source and lens distribu¬ 
tions. The result is 

p 00 poo 1 

<7 t(0))=l dzpi(z)AT,(R,z) dz p 3 (z)— -^-,(10) 

JO Jz Zjc(Z, Z ) 

noting that for narrow lens and source redshift distributions, 
Pi(z) = 5 d(z — z{) and p s (z) = 5 d(z — z„), Equation [TO] 
reduces to Equation [5] 

Equation m clarifies that for a broad lens redshift dis¬ 
tribution it is not possible to solve for AE (R,zi) in closed 
form. However, for a narrow lens redshift distribution, aver¬ 
aging over a background source redshift distribution p s (z), 
we can recover the relation 


AE (R, zi ) 


<7tW> 



dz p s (z) 


1 

E c (zi, z') 


-1 


( 11 ) 


Section Tf.21 presents an analysis to justify that this approxi¬ 
mation is acceptable, given the source and lens redshift dis¬ 
tributions of our datasets. Our method of estimating AE(7?) 
from the data is described in Section roi 


2.3 Galaxy clustering and redshift-space 
distortion 

The cross-power spectrum of lens galaxies and underlying 
mass that appears in Equation [3] P gm , and the equivalent 
cross-correlation function in Equation]!)] £ g m, depend on the 
manner or “bias” with which the lens galaxies trace the mat¬ 
ter field. Although this bias is in general a stochastic, non¬ 
linear and scale-dependent function, it may be approximated 
on sufliciently-large scales as a linear mapping between the 
galaxy and matter overdensity, 5 g (x) = b5 m (x) (e.g., Scher- 
rer & Weinberg 1998). In this case ^ gm (r) = b£mm(r), in 
terms of the matter auto-correlation function ^ mm (r) which 
may be modelled from the cosmological parameters, where 
the unknown bias parameter b may be determined using the 
auto-correlation function of the lenses, £ gg (r) = 6 2 ^ mm (r). 

Since 3D measurements of galaxy clustering are per¬ 
formed in redshift-space, the apparent radial positions of 
the galaxies contain an additional correlated imprint from 
galaxy peculiar velocities, known as redshift-space distortion 
(RSD). In particular, the Fourier transform of the redshift- 
space galaxy overdensity field, J g , is given under certain as¬ 
sumptions (e.g. Percival & White 2009, Blake et al. 2011) 
by 

s s g (k,p) = s g (k)-p 2 e(k), ( 12 ) 

where 6(k) is the Fourier transform of the divergence of the 
peculiar velocity field in units of the co-moving Hubble ve¬ 
locity and p, is the cosine of the angle of the Fourier mode 
to the line-of-sight. 

Assuming that the velocity field is generated under 
linear perturbation theory then 9(k) = —f 8 m {k ), where 
/ is the growth rate of structure, expressible in terms of 
the growth factor D(a) at cosmic scale factor a as / = 
din D/d In a. The growth factor is defined in terms of the 
amplitude of a single perturbation as 5(a) = D(a ) 5(1). Un¬ 
der the assumption of linear galaxy bias we then obtain 
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the standard expression for the redshift-space galaxy power 
spectrum in linear theory (Kaiser 1987) 

= 6 2 (l+/3 M 2 ) 2 P mm (fe), (13) 

where we have introduced the redshift-space distortion pa¬ 
rameter /3 = f/b, which governs the amplitude of the mea¬ 
sured RSD. 

The anisotropic imprint of RSD in galaxy clustering al¬ 
lows the measurement of the gravitational growth rate and, 
consequently, powerful tests of gravitational physics. How¬ 
ever, it also introduces an extra amplitude factor in the rela¬ 
tion between £ 99 and £ mm , complicating inferences about the 
galaxy bias. In order to avoid this issue the real-space “pro¬ 
jected” correlation function w p (R), independent of RSD, can 
instead be constructed by integrating the 3D galaxy corre¬ 
lation function £ 99 (.R, n) along the line-of-sight: 

/ OO 

£ 99 (-R,n)dn. (14) 

-oo 

Our method of estimating w p (R) from the data is described 
in Section IQ In practice the limits of Equation [14] must be 
taken as large, finite values. 


2.4 Suppressing small-scale information 


Equation [ 6 ] demonstrates that the amplitude of A 'S(R) de¬ 
pends on the surface density of matter around galaxies 
across a range of smaller scales from zero to R. This is prob¬ 
lematic from the viewpoint of fitting cosmological models 
to the data since at small scales, within the halo virial ra¬ 
dius, the cross-correlation coefficient between the matter and 
galaxy fluctuations is a complex function which is difficult 
to predict from first principles (Baldauf et al. 2010, Man- 
delbaum et al. 2010). In order to remove this sensitivity to 
small-scale information these authors proposed a new statis¬ 
tic, the annular differential surface density (ADSD), denoted 
by T and defined by 


T grn (R, Ro) = 


A£(P) - A£(Po) 

2 f R 

= jjj J R'nR')dR' 


- E(P) + ffS(Po), 


(15) 


which does not contain information originating from scales 
R < Ro- The small-scale limit Ro is chosen to be large 
enough to reduce the main systematic effects discussed 
above, but smali enough to preserve a high signal-to-noise 
ratio in the measurements (aiso see the discussion in Man- 
delbaum et al. 2013). An alternative approach is to model 
the halo occupation statistics and marginalize over the free 
parameters (e.g., Cacciato et al. 2013). 

The corresponding quantity suppressing the small-scale 
contribution to the galaxy auto-correlations is 


99 (^; Ro) — p . 

2 _ 

R 2 


r R 2 

/ R'w p (R')dR'-w p {R) + ^w p (R 0 ) 

JRq K . 


(16) 


We assume Ro = 1.5 6 _1 Mpc as our fiducial value following 
Reyes et al. (2010), and demonstrate in Section[6]3]that our 
results are insensitive to this choice. 


2.5 Testing gravitational physics: the Eg statistic 

In general scalar theories of gravity, the perturbed FRW 
spacetime metric ds 2 may be expressed in terms of the New¬ 
tonian potential 4/ and curvature potential 4>: 

ds 2 = [1 + 2 4t(x, t)] c 2 dt 2 — a{t) 2 [1 — 2 &(x, t)] dx 2 . (17) 


Relativistic particles, such as photons experiencing gravi¬ 
tational lensing, collect equal contributions from these two 
potentials as they traverse spacetime, such that their equa¬ 
tions of motion (and hence the resulting lensing patterns) 
are determined by V 2 (4t + 4>). However, the motion of non- 
relativistic particles arising from the gravitational attraction 
of matter, which produces galaxy clustering and RSD, is 
sensitive only to the derivatives of the Newtonian potential 
V 2 'E (e.g., Jain & Zhang 2008). 

In standard General Relativity (GR), in the absence of 
anisotropic stress, = &(x, t) and both potentials are 

related to the matter overdensity via the Poisson equation 
V 2 4> = 47rGo 2 p7TJm. Therefore, by measuring if both the 
gravitational lensing of photons and galaxy peculiar veloc¬ 
ity respond in an identical manner to the matter overdensity 
traced by the lens galaxies in our datasets, we can perform 
a fundamental test of whether the relation between (Jt 4 - <[>) 
and 4> follows the GR expectation (assuming this perturba¬ 
tion approximation applies). 

Zhang et al. (2007) proposed that this test can be effi¬ 
ciently carried out by cross-correlating lens galaxies to both 
the surrounding velocity field using RSD and to the shear of 
background galaxies using galaxy-galaxy lensing. In partic¬ 
ular, Reyes et al. (2010) implemented this consistency test 
by constructing the “gravitational slip” statistic 


E g (R) 


1 T gm(R, Ro) 
H ?gg(R,Ro) ’ 


(18) 


which is independent of both the galaxy bias factor b and 
the underlying amplitude of matter clustering as, given that 
oc 1/6, T grn oc bug and X 99 oc b 2 as- The perturbed GR 
model prediction on large scales is then a scale-independent 
quantity Eg = fi m // (see Leonard et al. (2015) for a more 
detailed discussion of this approximation). We measure Eg 
and carry out this consistency test in Section [6.31 We note 
that a failure of this consistency check does not necessarily 
indicate evidence for gravitational physics beyond GR: other 
possible explanations would include a breakdown in validity 
of linear perturbation theory, or that the value of or 
curvature differs from that predicted by measurements of 
the Cosmic Microwave Background radiation. 


3 DATA 

We perform this test of gravitational physics by utiliz¬ 
ing the overlap of lensing measurements from two imaging 
surveys, the Canada-France-Hawaii Telescope Lensing Sur¬ 
vey (CFHTLenS; Heymans et al. 2012) and the Red Clus¬ 
ter Sequence Lensing Survey (RCSLenS, Hildebrandt et al. 
2015), with two spectroscopic-redshift large-scale structure 
surveys, the WiggleZ Dark Energy Survey (Drinkwater et 
al. 2010) and the Baryon Oscillation Spectroscopic Survey 
(BOSS, Eisenstein et al. 2011). Figure [l] displays the sky dis¬ 
tribution of the CFHTLenS, RCSLenS, WiggleZ and BOSS 
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BOSS DR1 0 NGP 



BOSS DR10 SGP 



Figure 1. (R.A., Dec.) distribution in the Northern Galactic Pole (NGP) and Southern Galactic Pole (SGP) of the datasets used in 
this analysis: the BOSS DR10 galaxy sample (black dots), the WiggleZ survey regions (red rectangles), the CFHTLenS fields (green 
rectangles) and the RCSLenS fields (blue rectangles). 


Table 1. Statistics for CFHTLenS and RCSLenS fields cross-correlated with WiggleZ and BOSS data. The effective (unmasked) areas 
are shown for the full source field, the pointings which contain a minimum of four filters such that photometric redshifts are available, and 
the pointings which overlap with the lens distribution (where two values are quoted for this latter area, they refer to WiggleZ/BOSS). 
The effective weighted source density as computed by Equation 1191 is listed as n e ff. The number of lenses contained in the overlapping 
areas of WGZLoZ, WGZHiZ, LOWZ and CMASS are displayed as -NwgzloZ? -^WGzrnz? A^lowz and Admass? respectively. Some fields 
have overlap with both WiggleZ and BOSS lens samples, which would result in potentially correlated measurements. In such cases the 
lens sample producing the lower signal-to-noise measurement, indicated with an asterisk, is excluded from the analysis. 


Field 

A e ff [all] 
[deg 2 ] 

A e ff [photo- z] 
[deg 2 ] 

A e ff [overlap] 
[deg 2 ] 

n eft 

[arcmin -2 ] 

WvGZLoZ 

■WvGZHiZ 

ATowz 

ADMASS 

CFHTLSW1 

56.0 

56.0 

48.7 

14.0 

- 

- 

1998 

3856 

CFHTLS W4 

17.6 

17.6 

17.6 

13.1 

- 

- 

832 

1711 

RCS 0047 

51.8 

37.2 

37.2 

5.4 

3029* 

4343* 

2273 

3735 

RCS 0133 

25.0 

13.3 

13.3 

4.5 

- 

- 

646 

1236 

RCS 0310 

60.9 

52.3 

52.3 

4.9 

4249 

6140 

- 

- 

RCS 1303 

12.9 

8.3 

3.7 

5.3 

- 

- 

128 

412 

RCS 1514 

60.0 

30.4 

5.8/30.4 

5.7 

152* 

448* 

1504 

3337 

RCS 1645 

22.8 

20.1 

20.1 

6.7 

- 

- 

1098 

2008 

RCS 2143 

65.5 

47.1 

47.1 

5.7 

7190* 

8752* 

2337 

4718 

RCS 2329 

36.0 

32.1 

19.0/32.1 

6.4 

494* 

882* 

1356 

3122 

RCS 2338 

57.3 

19.5 

19.5 

4.9 

1336 

2158 

- 

- 


datasets used in this analysis, and the surveys and source 
selection are briefly described in the sub-sections below. 

A total of 11 CFHTLenS and RCSLenS survey fields 
overlap with the WiggleZ and BOSS DR10 data, compris¬ 
ing a total area of 466 deg 2 (74 deg 2 for CFHTLenS and 392 
deg 2 for RCSLenS). Table H lists some statistics for these 
fields, including the total effective (unmasked) field area, 
the area for which the available imaging allows photometric 


redshifts to be derived, the subset of that area which over¬ 
laps with the lens distributions, the effective source density 
(defined below) and the number of lenses in each of the over¬ 
lapping spectroscopic surveys used in the analysis, where 
the BOSS data is split into the LOWZ and CMASS sam¬ 
ples (described below). The RCSLenS fields used for cross¬ 
correlation with (WiggleZ, BOSS) contain an overlapping 
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photo -2 area of (72, 184) deg 2 . The CFHTLenS fields used 
for cross-correlation with BOSS cover 74 deg 2 . 

Our datasets enable us to construct five distinct source- 
lens survey pairings: RCSLenS-WiggleZ, RCSLenS-LOWZ, 
RCSLenS-CMASS, CFHTLenS-LOWZ and CFHTLenS- 
CMASS. We split the WiggleZ lenses into two indepen¬ 
dent redshift bins, 0.15 < 2 < 0.43 (“WGZLoZ”) and 
0.43 < 2 < 0.7 (“WGZHiZ”), which coincide with the red- 
shift ranges of LOWZ and CMASS, respectively, producing 
a total of six possible pairings. In the analyses that follow 
we will often present separate measurements for these six 
cases, optimally combining the measurements in each indi¬ 
vidual field using inverse-variance weighting. 


3.1 CFHTLenS 

The CFHTLenfJ]] is a deep multi-colour survey optimized 
for weak lensing analyses, observed as part of the Canada- 
France-Hawaii Telescope (CFHT) Legacy Survey in five op¬ 
tical bands u* g'r'i'z' using the 1 deg 2 camera MegaCam. 
The data span four fields, two of which (W1 and W4) over¬ 
lap with the spectroscopic-redshift data used in this anal¬ 
ysis. The main lensing analysis is performed on %' band 
data, for which the observations have a 5er point-source 
limiting magnitude i! « 25.5. The imaging data are pro¬ 
cessed with THELI (Erben et al. 2013), galaxy ellipticities 
are measured by the Bayesian model-fitting software Zensfit 
(Miller et al. 2013), and photometric redshifts are derived 
from PSF-matched photometry (Hildebrandt et al. 2012) 
using the Bayesian photometric redshift code BPZ (Benitez 
2000). Additive and multiplicative shear calibration correc¬ 
tions have been derived as a function of galaxy size and 
signal-to-noise (Heymans et al. 2012, Miller et al. 2013). 
The survey pointings have been subjected to a stringent 
cosmology-independent systematic-error analysis (Heymans 
et al. 2012), as a result of which a subset of around 25% 
of the pointings have been flagged as possessing potential 
systematic errors. Given that galaxy-galaxy lensing is much 
less sensitive than cosmic shear to such systematics due to 
the circular averaging over lens-source pairs, as thoroughly 
investigated for CFHTLenS by Velander et al. (2014), we 
do not remove these pointings from our analysis. We explic¬ 
itly checked that any difference in the galaxy-galaxy lensing 
statistics between these choices was within the range of sta¬ 
tistical fluctuations. 

The Zensfit pipeline returns measured ellipticity com¬ 
ponents (ei,e 2 ) for each source, together with an approx¬ 
imately optimal weight w 3 , a combination of the variance 
in the intrinsic ellipticity and the measurement error due to 
photon noise (Miller et al. 2013). We note that the sign of 
e 2 listed in the source catalogues must be reversed in our 
cross-correlation pipeline, because the positive ^-direction 
of pixel coordinates lies in the negative R.A. direction. We 
cut the catalogue to sources with weights w s > 0 which lie 
in unmasked areas of data. No magnitude cuts are applied, 
although fainter galaxies are downweighted by lower values 
of w 3 . We also do not apply a star-galaxy separation cut, 
since stars are already assigned negligible weights by Zensfit. 


The effective source density for lensing analyses (following 
Heymans et al. 2012) is defined by 


^eff 


1 

A e ff 


( E ,<) 2 

Ei «) 2 ’ 


(19) 


where A e fr is the effective (unmasked) area. The values de¬ 
rived for (Wl, W4) are n e fj = (14.0, 13.1) arcmin~ 2 . (We 
note that Chang et al. 2013 provide another possible defini¬ 
tion of n e ff.) 

The BPZ photometric-redshift pipeline returns a full 
redshift probability distribution pbpz( 2 :) for each source, 
binned in intervals of A 2 = 0.05 in the range 0 < 2 < 3.5, 
and we make use of this full information when computing 
the lensing signal. For purposes of binning sources in photo- 
2 slices, we use the maximum-likelihood redshift value zb of 
these distributions. The source redshift distribution p 3 (z), 
which is required for cosmological modelling and for build¬ 
ing the mock catalogues below, is constructed by stacking 
the pbpz ( 2 ) distributions for each source, weighting by the 
Zensfit weights, 


Ps ( 2 ) oc y^ y W 3 PBPZ ,i(z) 
i 


( 20 ) 


We do not apply any photo -2 cuts in our fiducial analysis, 
and demonstrate that our results are insensitive to such cuts. 

We note that Miyatake et al. (2015) and More et al. 
(2015) recently presented related galaxy-galaxy lensing mea¬ 
surements for CFHTLenS-CMASS, with the aim of under¬ 
standing the properties of the dark matter haloes traced by 
the lenses, and performing joint cosmological and halo oc¬ 
cupation model fits. 


3.2 RCSLenS 

The 2nd Red Sequence Cluster Survey (RCS2, Gilbank et 
al. 2011) is a nearly 800 deg 2 imaging survey in three optical 
bands g'r'z' carried out with the CFHT. The primary survey 
area is divided into 13 well-separated patches on the sky, 
each with an area ranging from 20 to 100 deg 2 . Nine of 
these fields (with properties listed in Table [l]) overlap with 
the spectroscopic-redshift data used in this analysis. The 
main lensing analysis is performed using the r' band data, 
with limiting magnitude r' « 24.3. Around two-thirds of the 
RCS2 area has also been imaged in the i' band. 

The RCS2 team have presented a series of investiga¬ 
tions of galaxy-galaxy lensing using these data, probing the 
occupation and shapes of dark matter halos (van Uitert et 
al. 2011, 2012, Cacciato et al. 2014) and the connection to 
galaxy luminosity, stellar mass and velocity dispersion (van 
Uitert et al. 2013, 2015). We focus instead on cosmological 
applications of this dataset. 

The RCSLenS datase10 is a lensing re-analysis of the 
RCS2 imaging data constructed by applying the same pro¬ 
cessing pipeline for shape measurement and photometric- 
redshift estimation as developed for CFHTLenS (Hilde¬ 
brandt et al. 2015). The effective RCSLenS source density 
derived using Equation 1191 n e fi ~ 5.5 arcmin -2 , is lower 
than the corresponding value for CFHTLenS due to the 


t http://www.cfhtlens.org 


1 http://www.rcslens.org 
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shallower imaging, but the data cover a significantly wider 
area. BPZ photometric redsliifts are only derived for point¬ 
ings containing 4 optical bands g'r'i'z' , which correspond to 
about two-thirds of the survey area (with statistics for indi¬ 
vidual fields listed in Tabic [lj. We note that the absence of 
u* band imaging in RCSLenS causes slightly larger photo-z 
errors and a greater outlier fraction at low redshifts, owing 
to more serious colour-redshift degeneracies. 

The RCSLenS shape catalogues have been subject to a 
“blinding” process to prevent confirmation bias in analysis 
(Hildebrandt et al. 2015, see also Kuijken et al. 2015 for a 
full description of the same process). In brief, the catalogues 
contain four sets of ellipticity data, one of which is the true 
data, whilst the other three sets have been manipulated by 
a trusted external colleague by different amounts unknown 
to the science team, sufficient to prevent confirmation bias 
but not so extreme as to render those data suspiciously dis¬ 
crepant. All scientific analyses are repeated for each of the 
4 blindings. Once the results are ready for publication, the 
blinding key is provided by the external, after which the sci¬ 
ence team can make no further undocumented modifications 
to analysis procedures. 

3.3 WiggleZ Dark Energy Survey 

The WiggleZ Dark Energy Survey (Drinkwater et al. 2010) 
is a large-scale galaxy redshift survey of bright emission-line 
galaxies over the redshift range z < 1 with median redshift 
z « 0.6, which was carried out at the Anglo-Australian Tele¬ 
scope at Siding Spring, Australia, between August 2006 and 
January 2011. In total, of order 200,000 redshifts were ob¬ 
tained, covering of order 1000 deg 2 of equatorial sky divided 
into seven well-separated regions, labelled as (00, 01, 03, 09, 
11, 15, 22) by their location on the sky in hours of R.A. 

WiggleZ galaxies were selected for observation using 
colour and magnitude cuts from a combination of optical 
and UV imaging. The optical imaging employed was from 
the SDSS in the Northern Galactic Pole and from the RCS2 
survey in the Southern Galactic Pole, spanning a total of 6 
RCS2 fields (0047, 0310, 1514, 2143, 2329, 2338). The Wig¬ 
gleZ dataset is therefore very well-suited for providing lenses 
for cross-correlation with RCSLenS sources. The 0047, 1514, 
2143 and 2329 fields are also covered by BOSS, producing 
higher signal-to-noise galaxy-galaxy lensing measurements 
in these cases, and therefore to avoid the complication of 
the covariance we focus on analyzing the 2 remaining fields, 
0310 and 2338. When measuring the galaxy-galaxy lensing 
signal, we also cut the lens and source distributions to the 
subset of MegaCam pointings containing both lenses and 
sources (the “overlap area”) listed in Table [T] After these 
cuts, we used a total of 5,585 WGZLoZ lenses and 8,298 
WGZHiZ lenses in our analysis. The WiggleZ lens selection 
function within each region was determined using the meth¬ 
ods described by Blake et al. (2010). The average galaxy 
bias factor of the WiggleZ lenses is b ss 1. 

Figure[2]plots the number density distribution with red¬ 
shift of the WiggleZ lens samples in the different survey re¬ 
gions. The typical number density is n « 3 x IIP 4 h 3 Mpc -3 
at the effective redshift z « 0.6, with a tail to lower number 
densities at high redshift (1 x 10 -4 h 3 Mpc -3 at z « 0.9) and 
higher number densities at low redshift (6 x 10 -4 h 3 Mpc~ 3 
at z ~ 0.1). Some variation amongst WiggleZ regions is ev- 


o 



z 


Figure 2. Number density distribution with redshift of lenses 
in each WiggleZ survey region and in the BOSS LOWZ/CMASS 
samples. The redshift distribution differs between WiggleZ re¬ 
gions because of varying colour/magnitude selection and com¬ 
pleteness of spectroscopic follow-up. 

ident, due to differences in the colour/magnitude selection 
and completeness of the spectroscopic follow-up (for full de¬ 
tails, see Drinkwater et al. 2010). 

3.4 BOSS 

BOSS is a spectroscopic follow-up survey of the SDSS III 
imaging survey (Eisenstein et al. 2011), which has obtained 
redshifts for over a million galaxies covering 10,000 deg 2 . All 
observations have been carried out at the Sloan Telescope 
located at Apache Point Observatory in New Mexico. BOSS 
uses colour and magnitude cuts to select two classes of galax¬ 
ies to be targetted for spectroscopy, the “LOWZ” sample 
which contains red galaxies at z < 0.43 and the “CMASS” 
sample which is designed to be approximately stellar-mass 
limited for z > 0.43. We use the data catalogues provided 
by the SDSS 10th Data Release (DR10); full details of these 
catalogues are given by Eisenstein et al. (2011), Dawson et 
al. (2013) and Anderson et al. (2014). 

Following the practice of the BOSS science papers we 
cut the LOWZ sample to 0.15 < z < 0.43 and the CMASS 
sample to 0.43 < z < 0.7, in order to avoid redshift overlap 
and create homogeneous galaxy samples. The redshift dis¬ 
tributions of these samples are shown in Figure [2] and are 
comparable to the WiggleZ dataset. The galaxy bias factors 
of the LOWZ and CMASS samples are b « 1.6 and b ~ 1.9, 
respectively (Chuang et al. 2014, Sanchez et al. 2014), hence 
these objects are significantly more clustered than WiggleZ 
galaxies. 

In order to correct for the effects of redshift failures, fi¬ 
bre collisions and other known systematics affecting the an¬ 
gular completeness (correlated with the stellar density and 
seeing), BOSS galaxies are assigned completeness weights 
(as specified in Equation 18 of Anderson et al. 2014). We 
included these weights in our determination of the auto¬ 
correlation function of the galaxies, but they do not affect 
the galaxy-galaxy lensing measurements because they are 
uncorrelated with the background shapes. 

As noted in Table [TJ the BOSS sample overlaps 2 
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Figure 3. Measurements of 7 1 ( 6 ) (left column) and 7 X ( 8 ) (right column) for the cross-correlation of different combinations of source-lens 
datasets. We show results for both the data (black solid circles) and mock mean (red open circles), with the errors based on measurements 
for a set of 374 mock catalogues. The overplotted model is the expectation for the mock mean based on the input cosmology of the 
simulations, x 2 statistics are quoted for the 7 X measurements with respect to a model of zero, with number of degrees of freedom 
dof = 20 . 


CFHTLenS fields (W1 and W4) and 7 RCSLenS fields (0047, 
0133, 1303, 1514, 1645, 2143, 2329). We utilized a total of 
12,172 LOWZ lenses and 24,135 CMASS lenses in our anal¬ 
ysis. 


4 GALAXY-GALAXY LENSING 
MEASUREMENTS 

In this Section we describe our estimation of the stacked 
tangential shear around the lenses as a function of angular 


separation and the differential surface density around the 
lenses as a function of the projected physical separation. In 
the latter case we use the full source photometric-redshift 
probability distribution when computing the lensing signal. 

4.1 Average tangential shear 7 t($) 

The ellipticity components (ei : i,e 2 ,i) of source i relative to 
the positive ir-axis are defined by 

(ei,e 2 ) = (^—ry) ( cos2 Vhsin2?/>), (21) 
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where tp is the position angle of a galaxy with axial ra¬ 
tio r measured anti-clockwise from the positive rr-axis. 
These ellipticity components can be rotated to new values 
[et(i, j), e x (*, j)] relative to a line connecting source i and 
lens j by 

= —ei : iCos2cj)(i,j) -e 2 ,jsin 2 <j>{i,j) ( 22 ) 

and 


4.2 Differential surface density AE (R) 

In Section we demonstrated that for a relatively nar¬ 
row lens distribution, AE(R, zi) could be estimated from 
the tangential galaxy shear via Equation 1111 which we re¬ 
write here as 

AE(J?,^) = 7tW/(sr T ), (26) 


ex(i,j) = ei,isin2</>(*,j) - e 2 ,icos20(i, j), (23) 


where (p(i,j ) is the angle of the line connecting source i 
and lens j to the positive x-axis (in the range —90° < <j> < 
90°). Our galaxy-galaxy lensing estimators for the tangential 
shear yt and cross shear -y x are then: 


7 1 { 0 ) = 
and 


7x(<?) = 


Esourcesi Elenses j W “ W ‘i ^ j) 9 (U) 
EsourcesiElen Bes j j) 

Esources £ Elen.es j W i W i ^ (b j) j) 

Esources* Elenses, W i W \) 6 (b j) 


(24) 


(25) 


where wf is the fensfit weight of source i, w l j is the weight 
of lens galaxy j which we set to the completeness weights 
for BOSS galaxies and w l = 1 for WiggleZ galaxies (since 
in that case the completeness is incorporated in the random 
catalogues), O (i,j) is equal to 1 if the angular separation 
between source i and lens j lies in bin 8 and equal to 0 
otherwise, and the sums are taken over unique pairs. We 
performed measurements in 20 logarithmically-spaced an¬ 
gular bins from 8 = 0.01° to 1°. Multiplicative shear bias 
corrections are computed and applied as discussed in Section 
EH The cross shear should be consistent with zero signal 
for weak gravitational lensing but potentially non-zero for 
distortions due to systematic errors. 

Figure |3] shows the resulting measurements of the tan¬ 
gential and cross shear for the six different combinations of 
source-lens datasets. The measurements are performed in 
the individual fields listed in Tabic [T] and then combined 
together using inverse-variance weighting. The errors in the 
measurements are obtained by applying the same analysis 
pipeline to a set of 374 mock catalogues, which are con¬ 
structed by the process described in Section[5] and appropri¬ 
ately scaling the standard deviation of the mock statistics. 
In Figure [3] we plot measurements for both the data (black 
solid circles) and the mean of the mock catalogues (red open 
circles), where the solid line is the model prediction for the 
mocks based on the input cosmology of the simulations. In 
the absence of systematic effects, the cross shear should be 
consistent with zero. Values of the chi-squared statistic are 
shown for each source-lens combination, determined for a 
model 7 x = 0 using the covariance matrix built from the 
mock catalogues. Our measurements are consistent with this 
expectation. 

At small scales, for values of 8 corresponding to pro¬ 
jected physical separations R < 2 h~ 1 Mpc, the amplitude 
of galaxy-galaxy lensing in the mock falls below that mea¬ 
sured in the data. This arises because of halo occupation 
effects on small scales particularly affecting Luminous Red 
Galaxies; it presents no concern for our analysis, in which 
physical effects originating at small scales R < Ro are ex¬ 
plicitly suppressed. 


where 


S C= dZsPs{Zs) E c (zi,z„) 

• z l 

AttG (1 + zi) x(zi) 


dz s Ps(z s ) 


1 - 


X(zi) 

X(zs) 


(27) 


We implemented this calculation for each lens-source pair 
by using the source photometric-redshift probability distri¬ 
bution, p s (z) = pbpz(z). Therefore we estimate AE using 
only that lensing data for which photometric-redshift infor¬ 
mation is available, as listed in Table [I] 

For our determinations of E7 1 using Equation [27] and 
conversions of source-lens angular separations 8 to projected 
separations R = 8x(zi), we adopted a fiducial flat ACDM 
cosmological model with matter density fl m = 0.27. Our 
motivation for this choice, which is in better agreement with 
the fits to the Cosmic Microwave Background fluctuations 
using WMAP satellite data (Komatsu et al. 2011) rather 
than the later Planck satellite data (Planck collaboration 
2015b), is to ensure consistency with the fiducial cosmologi¬ 
cal model adopted for the RSD analyses of the WiggleZ and 
BOSS data (Blake et al. 2011, Sanchez et al. 2014), which 
would be subject to Alcock-Paczynski distortion in different 
models. Adopting the higher value of preferred by the 
Planck analysis would not produce a significant change in 
our measurements compared to the statistical errors. 

When averaging the estimates of Equation [7S] over the 
lens-source pairs, each source must now be inverse-variance 

weighted as w s ^EJ 1 ^ and hence our final estimator is 


AE (R) = 


^—/sources i i le 


j ^7 Wj E c jj et{i,j) 0(1 , j j 


^sources i / 


lenses j W i W j 


(28) 


O (i,j) 


(see also Miyatake et al. 2015). This estimator corrects for 
the dilution in lensing signal caused by the non-zero proba¬ 
bility that a source is situated in front of the lens, z 3 < zi. 

We note that an approximate version of this estimator 
has been used in previous studies where only the maximum- 
likelihood source redshift z 3 is available, rather than the 
full probability distribution p s (z s ). In this approximation 
the lens-source pairs are restricted to those cases for which 
z., > zi and the differential surface density for each pair is 
estimated as 


AE (R, zi) 7 1 {8) E c (zi, z 3 ). 


(29) 


The inverse-variance weight of each source-lens pair is now 
w s E “ 2 such that the approximate estimator reads 

E -Ei . wl ErJ,. e t(i, j) Q(i,j) 

* v-i/7->\ /sources z ^—/lenses j z J c,xj L \ ’ J / \ ’ ^ / / or »\ 

ZAZj(ixJ ^ ^ " / v ^-2 r-\t ■ ~ ■ V’’ ) 


^—/sources i ^—/lenses j 


. wf w l ; E ?■ Q(i,j) 
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Figure 4. AE(i?) measured for the cross-correlation of different combinations of source-lens datasets. We show results for the data and 
the mock mean. The overplotted model is included only for comparison with the mock mean. 
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Figure 5. A comparison of the measurements of AS (R) for different combinations of source-lens datasets using the approximate estimator 
of Equation |30| based only on the best source photometric redshift zb , plotted as the solid (black) circles, and the unbiased estimator of 
Equation 1 28 1 based on the full redshift probability distribution PbpzW °f each source, plotted as the open (magenta) circles. For the 
purposes of this plot the errors were determined by jack-knife re-sampling. 
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- WGZLoZ—RCSLenS z eff =0.32 ; 

-LOWZ—RCSLenS z eff =0.32 

-LOWZ-CFHTLenS z eff =0.32 

WGZHiZ—RCSLenS z -=0.57 



Figure 6. Verification that the estimator of Equation 1281 suc¬ 
cessfully recovers the input value of AE (R) in mock catalogues. 
A toy photo-z model is applied to the mocks, in which redshifts 
are scattered in accordance with a Gaussian distribution of zero 
mean and standard deviation 0.2 x (1 + z). The differential sur¬ 
face density is then estimated using both Equation 1301 restricting 
the calculation to source-lens pairs for which z p hot > ^spec (red 
open circles), and Equation 1281 which uses the full information of 
the photo-z probability distribution (green open triangles). These 
determinations are compared to the measurements without appli¬ 
cation of photometric redshifts (black solid circles). In all cases 
the mock mean is plotted, where the error represents the standard 
deviation of the mocks. 


This relation will contain a bias as far as E c ^ 1/E^" 1 , al¬ 
though that bias can be corrected through comparison with 
spectroscopic catalogues (Nakajima et al. 2012). 

Figure 3] shows the measurement of AE (R) using Equa¬ 
tion m in 20 logarithmically-spaced bins in R from 0.5 to 
50 f* _1 Mpc, for the six different combinations of source-lens 
datasets, displaying results for both the data and the mock 
mean. As above, the errors in these measurements are ob¬ 
tained using the 374 mock catalogues introduced in Section 
0 Our choice of this binning scheme allows us to search for 
scale-dependence in the measurement whilst retaining good 
signal-to-noise in each bin. 

We now consider in more detail the difference in re¬ 
sults produced by the unbiased estimator of Eauation l28l and 
the approximate estimator of Equation l30l First, Figure [5] 
compares measurements of AE(.R) using the two estimators. 
We conclude that the approximate estimator contains a sig¬ 
nificant systematic bias, especially when using photometric 
redshifts with precision characteristic of the RCSLenS data. 
Secondly, Figure[S]presents a verification of our successful re¬ 
covery of AE (R) in the presence of photometric-redshift er¬ 
rors, using the mock catalogues described in Section[5] A toy 
photo -2 model is applied to the mock catalogues, in which 
redshifts are scattered in accordance with a Gaussian distri¬ 
bution of zero mean and standard deviation 0.2 x (1 + 2 ), 
and the differential surface density is then estimated using 
both Equation m and Equation m The measurements in 
which the full source redshift probability distributions are 
used agree well with an analysis in which the photo -2 scat¬ 
ter is not applied. 

Finally, Figure [7] assesses how accurately the narrow 
lens approximation described in Section m holds for the 


Figure 7. The systematic error in determination of AE that 
results in the application of the narrow-lens approximation of 
Equation ! 1 H at the effective redshift of the RSD measurements of 
each lens sample. 


specific lens and source redshift distributions of our sam¬ 
ples. The amplitude of the systematic error is calculated by 
comparing the value of AE (R) that would be inferred at the 
effective redshifts at which RSD is measured for the differ¬ 
ent lens samples by substituting Equation ! lOl in Equatioii llil 
with the fiducial value of AE (R) evaluated using Equation^ 
We find that the resulting amplitude of the systematic error 
in AE (and consequently Eg) is about 5%, which is signif¬ 
icantly smaller than the statistical error in Eg of around 
20%. We conclude that the narrow-lens approximation is 
acceptable for our analysis. 

We note (without further calculation) an additional is¬ 
sue for estimating Eg that could occur for wide lens red¬ 
shift slices. The weighting applied to each source-lens pair 
when estimating AE (R), EJ 1 , contains a dependence on 
the lens redshift which may induce a systematic difference 
in the relative weighting of lens galaxies to that used in 
the measurement of their projected auto-correlation func¬ 
tion w p (R), which would become important if the bias of the 
lenses evolved across the redshift slice. However, the bias of 
both WiggleZ and BOSS lenses is a fairly flat function of 
redshift across these ranges (Blake et al. 2009, Anderson et 
al. 2014). 


4.3 Calibration and random-catalogue corrections 

Bias in ellipticity measurements may be described by a linear 
combination of a multiplicative error m and an additive error 
c such that 

e° bs = (1 + mi) ef ue + Ci + Noise, (31) 

where * = 1,2 denotes the two ellipticity components. 

Velander et al. (2014) show that the additive shear bias 
has a negligible effect for galaxy-galaxy lensing measure¬ 
ments with CFHTLenS data. We follow the standard prac¬ 
tice of performing this calibration by subtracting the cor¬ 
relation with a random lens catalogue as described below, 
which represents an empirical calculation of any additive 
shear residual. 
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Figure 8. The cross-correlation between shapes and random lenses, AH r;IT1( ]( R). measured for the different combinations of source and 
lens samples. The results are averaged over 40 random catalogues. A jack-knife error is plotted, noting that this error is only accurate 
for small scales and is a significant over-estimate for the largest scales. This systematic correction is subtracted from the data, which is 
standard practice for calibrating galaxy-galaxy lensing measurements (see for example Mandelbaum et al. 2005). 


The multiplicative shear bias correction m defined in 
Equation [BT] has been modelled for lensfi t processing of 
MegaCam data using a shear recovery test based on galaxy 
image simulations as a function of source signal-to-noise SN 
and size in pixels r: 


m{SN ’ V) = lo gl0 %N) 6XP ( “ r SN a) ’ (32) 

where a = 0.057 and /3 = —0.37 (Miller et al. 2013). The 
same correction applies to both ei and e2, and is greater 
for sources with smaller signal-to-noise ratios and sizes. We 
propagated this correction into the measurement of 7 t {9) by 
evaluating 


K(9) 


Esourcesi Elensesj W t W i Q(m) 
Esourcesi Elensesj W i W \ 6 (b j) 


such that the corrected measurement is given by 


corrected 

it 


( 0 ) 


uncorrected 

7 1 


( 0 ) 


1 + K{6) 


(33) 


(34) 


We note that this multiplicative shear bias correction must 
be applied in a global fashion, not on an individual source 
basis, due to its correlation with the measured source el- 
lipticity values owing to its dependence on SN and s in 
Equation [32] 

The analogous formula for the correction to be applied 
to AE (R) is 


K(R) 


Esources i Elenses j W i W ‘i 9(«.j) 

Esourcesi Elensesj W i W J ^ c,ij ®(*) j) 


(35) 


Equation[35]assumes that we are using the approximate esti¬ 
mator of Eauation l30l when instead employing the unbiased 
estimator of Equation 1281 we should replace Eby EJ 1 ^. 
Figure IA2I in Appendix IA2I displays the corrections mea¬ 
sured for the different source and lens samples; we find that 
K « —0.06 independent of scale, such that the amplitude of 
the galaxy-galaxy lensing signal must be boosted by ~ 6%. 

The shear estimated by these methods might still be 
contaminated by two further effects: (1) large-scale resid¬ 
ual shape-measurement systematics, such as an imperfectly- 
modelled optical distortion across the camera field; (2) the 
physical association of source and lens galaxies, for exam¬ 
ple sharing the same dark matter halo, diluting the cross¬ 
correlation signal since these sources will not be lensed and 
yet may have been scattered to higher redshifts by photo-;? 
error. The importance of both of these effects can be deter¬ 
mined using random lenses sampled from the same selection 
function as the data, where we averaged over 40 random 
catalogues in these analyses (checking that our final results 
were insensitive to the number of catalogues and the uncer¬ 
tainty in the correction). 

First, we re-ran the shear measurement replacing the 
data lenses with the random lens catalogues. The result for 
the differential projected surface density, which we denote 
AE ran d(i?), is displayed in Figure [8] as the combined signals 
for the different combinations of source and lens samples. 
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In the absence of residual coherent tangential distortions 
the signal should be consistent with zero; however, a sig¬ 
nificant correction is obtained at large scales for some of 
our samples, in particular CFHTLenS, and is always sub¬ 
tracted from our measurements of AE (R). We investigate 
the effect of physically-associated sources in Appendix IA51 
concluding that it is not significant for the scales we study. 
We compared our results to those presented by Miyatake et 
al. (2015). These authors find a more significant large-scale 
signal in AE ran d ( R ) for CMASS-CFHTLenS than we do, al¬ 
though we note that they analyze a different geometry of 
CMASS lenses. 


4.4 Summary of systematics tests 

We supported our science measurements with a series of 
tests for potential systematic errors. These results are de¬ 
scribed in detail in the Appendices, and we provide a brief 
summary in this Section. 

In Appendix IAll we present a series of systematics tests 
based on re-measuring 7 1 (0) following various manipulations 
of the shear catalogue: 

• Rotation of the sources by 45° (i.e., ei, ne w = e 2 ,oia and 

62 ,new — 6l,old)- 

• Randomizing the shear values amongst the source cat¬ 
alogue. 

• Replacing the lens catalogue by a random catalogue. 

The results of these tests are consistent with yt = 0, as 
expected. 

Systematic errors in the BPZ photometric redshift dis¬ 
tributions would imprint errors in the determination of 
AE (R) (given that these full probability distributions are 
used in the measurement, as described above). In order to 
search for such effects we performed a “scaling test” of the 
galaxy-galaxy lensing signal measured for the same set of 
foreground lenses using background sources in a series of 
different photo- 2 : slices spanning the range 0 < zb < 1 - 6 . 
The results are detailed in Appendix IA3I We find that a 
consistent lens singular isothermal sphere velocity disper¬ 
sion is obtained for all source photo -z slices with the possi¬ 
ble exception of 0 < zb < 0 . 2 , which comprises a negligible 
number of galaxies. As a final test, we repeated the Eg mea¬ 
surement using sources in different zb ranges; these results 
are discussed in Section RUil and demonstrate that the mea¬ 
surement is robust to these choices. 

Another potential source of systematic error is intrinsic 
alignment of the source population with respect to the fore¬ 
ground lenses (Blazek et al. 2012), which is expected to pref¬ 
erentially diminish the average tangential shear of red source 
galaxies. In Appendix IAll we carry out shear measurements 
for red and blue sub-samples to verify that the amplitude 
of these effects are below the level that would impact our 
results; we will demonstrate this explicitly in Section [S3] by 
repeating the Eg measurement for a blue sub-sample. 

I 11 addition to the systematic tests carried out in the 
paper, we refer the reader to Kuijken et al. (2015) who apply 
the same data analysis software used in this analysis to data 
from the Kilo Degree Survey. 



8 [deg] 


Figure 9. Comparison of the mean of various clustering and lens¬ 
ing statistics measured for the mock catalogues, with the expec¬ 
tations based on the input cosmological model of the simulation 
and the source and lens redshift distributions and bias factors. Re¬ 
sults are shown for the lens and source auto-correlation functions 
“’lens and resource, the cross-correlation function w SOUICe _ i ena , the 
cosmic shear statistics 5 + and §_, and the galaxy-galaxy lensing 
signal 7 1 - The error bars correspond to the standard deviation 
across the mock catalogues (which do not contain shape noise). 


5 SIMULATIONS 
5.1 Mock catalogues 

In order to test our analysis pipeline and estimate the co- 
variance matrix of the measurements, we created mock cat¬ 
alogues of each source and lens field. These catalogues were 
built from a set of 374 N-body simulations created us¬ 
ing methods similar to those described by Harnois-Deraps, 
Vafaei & van Waerbeke (2012). In brief, the N-body simu¬ 
lations are produced by the CUBEP 3 M code with a transfer 
function obtained from CAMB (Lewis, Challinor & Lasenby 
2000 ) using the following cosmological parameter set: matter 
density fl m = 0.2905, baryon density fib = 0.0473, Hubble 
parameter h = 0.6898, spectral index n s = 0.969 and nor¬ 
malization as = 0.826. The box-size of the simulations is 
L = 505 h~ x Mpc; this is significantly larger than the simu¬ 
lation set used for modelling the earlier CFHTLenS measure¬ 
ments [L = (147, 231) /i -1 Mpc] such that the new simula¬ 
tion set is much less affected by suppression of the large-scale 
variance by finite box size. For each simulation the density 
field is output at 18 snapshots, and a survey cone spanning 
60 deg 2 is constructed by pasting together these snapshots, 
where the division between adjacent snapshots is taken as 
the mid-point of the cosmic distances corresponding to the 
output redshifts. The two-component shear fields are also 
computed at each output redshift, by ray-tracing through 
the survey cone using Born’s approximation. More details 
are provided by Harnois-Deraps & van Waerbeke (2015). 

For every different survey region and lens-source sur¬ 
vey pairing, we converted each simulation line-of-sight into 
a mock catalogue by generating: 

• A source distribution with the surface density values 
given in Table |T] (adjusted for each region), using the RC¬ 
SLenS or CFHTLenS source redshift distribution as appro¬ 
priate, and Monte-Carlo sampling sources from the density 
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field with bias b source = 1. The source shear components 
are assigned by linearly interpolating the shear fields at the 
source redshift from the enclosing lens planes. 

• A lens distribution with the numbers given in Ta¬ 
ble [T] (adjusted for each region), using the average Wig- 
gleZ, LOWZ or CM ASS redshift distribution as appropri¬ 
ate, sampled from the density fields with bias factors &i ens = 
(1.0,1.6,1.9) for WiggleZ, LOWZ and CMASS, matching 
the bias measurements of Blake et al. 2010, Chuang et al. 
2014 and Sanchez et al. 2014, respectively. 

Sources and lenses are produced with a continuous distri¬ 
bution in redshift by interpolating across the finite redshift 
width of each simulation snapshot; at this stage the angular 
selection function is uniform across the 60 deg 2 cone. For a 
given bias factor b, the galaxy density field p g is related to 
the mass density field p m as 


We note that a bias model with b > 1 cannot be applied 
self-consistently on all scales because then p g contains neg¬ 
ative regions. We avoid these regions by imposing the con¬ 
dition p g = max(p s , 0 ), which reduces the effective value of 
the large-scale bias by a few per cent (we settled on this 
approach after experimenting with various approaches for 
ameliorating this effect such as smoothing the density field 
before applying the sub-sampling). 

Figure [9] tests the resulting mock catalogues by mea¬ 
suring the average tangential shear 7 1 , cosmic shear signal 
(£+,£-) and source-lens auto- and cross- angular correlation 
functions w(9) for the RCSLenS-WiggleZ mock catalogues 
(prior to the application of shape noise). The lines are the 
model predictions assuming the input cosmological param¬ 
eters of the simulation, and source and lens redshift dis¬ 
tributions and bias factors. The matter power spectrum is 
generated using the CAMB transfer function with non-linear 
contribution from HALOFIT (Smith et al. 2003). The mock 
mean closely follows the model predictions in each case. 

Shape noise is applied to the source catalogues using 
the following method: 

• A complex noise n = m +ri 2 i is formed for each source, 
where ru and 17,2 are drawn from Gaussian distributions with 
standard deviation <r e . 

• A complex shear 7 = 7 ! + 72 i is formed from the shear 
components ( 71 , 72 ) obtained from the ray-traced shear field 
at the source position, linearly interpolating between the 
values at adjacent snapshot redshifts. 

• A complex noisy shear is then determined as e = (7 + 
n)/(l + « 7 *) (Seitz & Schneider 1997). The components of 
the observed shear (ei, 62 ) are then found as e = e\ + e 2 i- 

We used the value a e = 0.29 for our mocks, which is repre¬ 
sentative of the two imaging surveys (Heymans et al. 2013, 
Hildebrandt et al. 2015) albeit slightly on the conservative 
side. The result of this procedure is 374 mock catalogues, 
matching the global properties of the source and lens sam¬ 
ples with a uniform angular selection function across 60 
deg 2 . 

We also constructed mock catalogues for each survey 
region including the full survey mask of sources and lenses, 
implemented by stitching together multiple simulation lines- 


of-sight and sub-sampling the result to match the survey se¬ 
lection functions. 149 independent mock catalogues for each 
survey region could be generated from the 374 simulations. 
This masked simulation set allows the importance of the 
survey selection function in the measurement error to be 
determined, as described in the next sub-section. 

5.2 Determination of the covariance matrix 

We compared several techniques for obtaining errors in the 
measurement of 7 1 ( 9 ) and AE (R)\ 

• Data jack-knife errors, where 16 (4 x 4) jack-knife re¬ 
gions of typical dimension ~ 2 ° are used per lensing survey 
region, obtained by dividing the source distribution into sub¬ 
samples containing equal number of galaxies using constant 
R.A. and Dec. boundaries. 

• Simulated errors not including the survey mask. We 
used the simulations discussed above, which comprise 374 
lines-of-sight each covering 60 deg 2 . We measured the cross¬ 
correlations for each line-of-sight, and scaled the resulting 
scatter by \J (60 deg 2 )/A e ff where A e fj is the effective (un¬ 
masked) area of each source region listed in Tabic [l] 

• Simulated errors including the survey mask, imple¬ 
mented by stitching together multiple lines-of-sight for the 
simulations and sub-sampling the result to match the se¬ 
lection functions of the sources and the lenses, as discussed 
above. We measured the cross-correlations for each of the 
resulting 149 mock catalogues per region; the standard devi¬ 
ation of these measurements constitutes our error estimate. 

• Jack-knife errors applied to both types of mock cata¬ 
logues. For the mocks without the survey mask, we scaled 
the resulting scatter by (60 deg 2 )/A e ft for each region. 

Figure [in] compares these error determinations for the 
AE (R) statistic for the different combinations of source-lens 
datasets, combining results in the different survey regions. 
The interpretation of these results is that on small scales, the 
errors in the measurements are dominated by shape noise. 
At larger scales other effects become important, such as the 
same source galaxies contributing to many stacks around 
different lenses, and the “sample variance” contributed by 
the particular realization of large-scale structure within our 
lens samples. 

All the different error estimates agree well for scales R < 
5 h -1 Mpc for the low-redshift lenses (0.15 < z < 0.43) and 
R < 8 / 1 -1 Mpc for the high-redshift lenses (0.43 < z < 0.7), 
following a scaling of the projected separation correspond¬ 
ing to a given angular separation. The jack-knife method 
applied to the simulations produces an error estimate which 
matches on all scales that obtained from the actual data. 
The error determined from the area-scaled unmasked simu¬ 
lations and fully-masked simulations matches for all angles, 
suggesting that the survey mask is of secondary importance 
for determining the measurement errors compared to the 
total unmasked area. However, these simulation errors are 
somewhat lower than the jack-knife errors for larger separa¬ 
tions. We expect the jack-knife errors to become unreliable 
for large scales as the jack-knife regions become increasingly 
less independent. 

As a result of these tests, we constructed our measure¬ 
ment errors using the ensemble of Weal = 374 mock cat- 
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- 1 


(36) 
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Figure 10. Comparison of the errors in AS(R) determined by jack-knife re-sampling of the data (labelled Data JK), by using the 374 
simulation lines-of-sight and scaling by an effective area factor (comparing jack-knife re-sampling of the simulations labelled Sim+scale 
JK, and the scatter between the realizations labelled Sim+scale real), and by generating 149 simulations of each region including the 
full selection functions (comparing jack-knife re-sampling labelled Sim+mask JK, and the realization scatter labelled Sim+mask real). 
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Figure 11. Covariance matrices of the AE(it) measurements for different combinations of source-lens datasets, determined using the 
scatter across the 374 simulation lines-of-sight. The covariance matrix Cij is displayed as a correlation matrix Cij / \JCuCjj. 
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alogues without the survey mask. Figure [IT] illustrates the 
off-diagonal elements of the resulting covariance matrix of 
the AE (R) statistic for the different combinations of source- 
lens datasets, constructed as 


Cov[AE(i?i), AS(f?.j)] = —j- x 

-^real 

Y AE \Ri) AE k {Rj) - AE (ft) AE (Rj) , 

. k= 1 


(37) 


where AE fe is measured for the fcth mock catalogue, and 
AE(7?i) = AE fc (i?i)/iV r eai. Significant correlations 

between bins are evident for R > 10 h -1 Mpc. 

The inverse of these covariance matrices is used in the 
parameter fits described below. We correct the inverse co- 
variance for the bias due to its maximum-likelihood estima¬ 
tion (Hartlap et al. 2007) via multiplication by the factor 


Weal “ Win 2 
-Weal ~ 1 


(38) 


where Win is the number of data bins used in the fit. For 
our analyses, Win/Weal ~ 0.05. 


6 COSMOLOGICAL RESULTS 

6.1 Projected galaxy auto-correlation function 

w p (R) 

We measured the 2D galaxy correlation function £, gg (R, II) of 
each lens sample, binning galaxy pairs by projected separa¬ 
tion R and line-of-sight separation II. We hence determined 
the projected correlation functions 

w p (R) = 2 Y Zsv ni ) An i, (39) 

bins i 

where we summed over 10 logarithmically-spaced bins in II 
from 0.1 to 60 /i _1 Mpc. The measurements of w p (R) in 20 
logarithmically-spaced bins in R from 0.5 to 50 W 1 Mpc for 
the different lens samples are plotted in Figure ITUl including 
jack-knife errors. These projected correlation function mea¬ 
surements are used for simple parameter fits in Section IPl 
and to determine the Eg statistic in Section IQ 





Figure 12. The projected correlation function w p (R) for the dif¬ 
ferent lens data samples used in our analysis. Jack-knife errors 
are plotted, together with the best-fitting model using both the 
w p (R) and AE( R) measurements. 


>s 


_Q 

O 

_Q 

O 

ol 



CT 


8 


Figure 13. The posterior probability distributions of the as pa¬ 
rameter, after marginalizing over the bias factors of the lens galax¬ 
ies, for the different combinations of source-lens datasets. We also 
show the probability distribution of the fit to all datasets. 


6.2 Measurements of as 

As an initial consistency test of the best-fitting cosmological 
parameters in the ACDM model, we fitted the measurements 
of AE (R) and w p (R) for each source-lens combination vary¬ 
ing just as and the galaxy bias of each lens sample, fixing 
the other cosmological parameters at the values used to con¬ 
struct the N-body simulations listed in Section[5] Given that 
AE oc b as and w p oc b 2 as, the degeneracy between these 
normalization parameters is broken and they can be sepa¬ 
rately determined. Fits were performed using the full covari¬ 
ance matrix of AE (R) determined from the mock catalogues, 
and a diagonal error matrix for w p (R) using jack-knife er¬ 
rors. The use of the latter is not significant: the signal-to- 
noise of w p (R) is much higher than that of AE(7?), and the 
limiting factor for the final parameter error is the noise in 
AE, given the initial degeneracy between b and as- 

We calculated model predictions using Equation [9] for 


AE(7?) and Equation 1141 for w p (R) and fit to the measure¬ 
ments over the range of scales R > 5 /i -1 Mpc (noting that 
our results do not depend significantly on the choice of min¬ 
imum fitted scale). The best-fitting models to the data are 
overplotted in Figures [12] and 1141 The value of the chi- 
squared statistic for the best-fitting model is \ 2 = 128.3 
for 113 degrees of freedom, indicating a good fit. 

The combined measurement is as = 0.79 ± 0.07 consis¬ 
tent with the latest determinations from the Planck satel¬ 
lite (Planck collaboration 2015b). Combining the separate 
fits to determine the normalizations in two independent 
redshift bins, 0.15 < 2 < 0.43 (WGZLoZ, LOWZ) and 
0.43 < 2 < 0.7 (WGZHiZ, CMASS), we determine a$(z = 
0.32) = 0.75 ± 0.08 and as{z = 0.57) = 0.54 ± 0.07. 

In Figure [13] we display the posterior probability dis¬ 
tributions of as for each source-lens dataset, marginalized 
over the bias factors. The datasets offer roughly comparable 
constraining power, with CMASS-RCSLenS and WGZHiZ- 
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Figure 14. The annular differential surface density statistic for the galaxy-mass cross-correlation, Y grn (R, Rq ), measured for the different 
combinations of lens-source datasets assuming Ro = 1.5 1 Mpc. We also plot the best-fitting model for each cross-correlation using 

both the w p (R ) and AE (R) measurements. The errors are based on measurements for a set of 374 mock catalogues. The horizontal 
dotted line marks Y grn = 0. 




R [h -1 Mpc] R [h -1 Mpc] 


Figure 15. Ec{R) measurements in two independent redshift bins 0.15 < z < 0.43 and 0.43 < z < 0.7, after combining the results 
from the different cross-correlations. In the former case, the measurements of Reyes et al. (2010) are plotted as the open circles for 
comparison. The horizontal solid lines are the prediction of standard gravity, Eq = 11,,,//, for our fiducial model = 0.27. The 
horizontal dotted lines indicate the 1-a variation that would result given Af2 m = 0.02, which is indicative of both the WMAP and 
Planck error in determining this parameter. We note that the data points are correlated, with a covariance matrix displayed in Figure 

HU 


RCSLenS producing the most and least accurate determi¬ 
nations, respectively. 

As a cross-check of the methodology we performed the 
same fits to the AS (R) measurements from the mock cat¬ 
alogues for all the combinations of source-lens datasets, us¬ 
ing the full-survey realizations including masks. The aver¬ 
age parameter measurement across the realizations is as = 


0.80 ± 0.03 with average value of x 2 /dof = 50.5/47, com¬ 
pared to the input parameter value <T8 = 0.826. The slight 
offset of the fit to lower values than the input is due to the 
artificial reduction in the clustering amplitude of high-bias 
mocks constructed via Equation M as discussed in Section 
[5] For 6=1 mocks we recover the input cosmology within 
the statistical error in the mean. 
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Future work will perform a full cosmological parameter 
analysis of these lensing and clustering datasets, in combi¬ 
nation with the CMB. 

6.3 Measurement of gravitational slip Eg{R) 

In this Section we use the measured galaxy-galaxy lens¬ 
ing cross-correlations, in combination with the clustering 
strength and redshift-space distortion properties of the 
lenses, to compute the scale-dependent Eg statistic defined 
in Equation M and hence carry out a test of gravitational 
physics. We calculated Eg(R) for each different combination 
of source-lens datasets, and then combined the Eg mea¬ 
surements in two independent redshift bins, 0.15 < a < 
0.43 (WGZLoZ, LOWZ) and 0.43 < z < 0.7 (WGZHiZ, 
CMASS). 

Firstly we converted the measurements of AS (R) and 
w p (R) for each source-lens combination into the annular dif¬ 
ferential statistics T grn (R) and T gg {R) defined by Equa¬ 
tions and 1161 respectively. We determined the values of 
AE(flfl) and w p (Ro) via a power-law fit to the appropri¬ 
ate statistic, taking the fitting range as Ro/3 < R. < 3Rq. 
Following Reyes et al. (2010) we assumed a fiducial value 
R 0 = 1.5 h^ 1 Mpc for our analysis, although we also con¬ 
sider a range of other choices between 1 and 3 /r 1 Mpc. We 
propagated the errors in AE(f?o) and w p (Ro) into the er¬ 
rors in T gm.(R) and T gg (R), although this source of error is 
negligible except when R « Ro- We determined the integral 
appearing in the first term of Equation [16] using a spline 
fit to the measured w p (R) (noting that these details make a 
negligible contribution to error propagation). Measurements 
of T grn (R, Ro = 1.5h -1 Mpc) are displayed in Figure fill for 
the different source and lens combinations. 

We cross-checked the errors in the determination of 
T g m(R) by repeating the procedure for the ensemble of 
mock catalogues. The scatter in the measurements across 
the mocks agreed with the propagated value of the error 
within a few per cent. 

We then determined Eg(R) for the different samples, 
using the values of /3 quoted in the respective redshift-space 
distortion analyses of the lenses, as reproduced in Table [2] 
For the BOSS galaxies we used /Jlowz = 0.38 ± 0.11 and 
/Icmass = 0.36 ± 0.06 (Sanchez et al. 2014). We performed 
new RSD fits to the WiggleZ data for redshift ranges 0.15 < 
2 < 0.43 and 0.43 < z < 0.7 following the procedures of 
Blake et al. (2011), obtaining /3wgzloZ = 0.66 ± 0.10 and 
/SwGZHiz = 0.63 ± 0.08. We propagated the errors in T gm , 
T gg and f5 into the error in Eg- We then combined the 
different measurements of Eg(R) in the two redshift bins 
0.15 < 2 < 0.43 and 0.43 < z < 0.7, using inverse-variance 
weighting. 

The final results are displayed in Figure fl5l correspond¬ 
ing to a ~ 20% determination of Eg at each scale and red¬ 
shift. We overplot the GR model prediction at each redshift, 
Eg = 0.m/f{z) for our fiducial model = 0.27 (together 
with its 1 -a variation give the typical error in determination 
of this parameter from the CMB, Afl m = 0.02), and the 
measurements of Reyes et al. (2010), which correspond to 
the lower-redshift bin. The covariance matrix of these Eg 
measurements is shown in Figure [TB] 

Using the covariance matrix, we fit a model of con¬ 
stant Eg(R) to the data, comparing fits over the ranges 


R > 2R 0 = 3 h -1 Mpc and R > 10/C 1 Mpc. The resulting 
values of Eg for each individual source-lens dataset are listed 
in Table [5] The combined measurements in the two unique 
redshift bins 2 = (0.32,0.57) are (0.40±0.09, 0.31 ±0.06) for 
R > 3h -1 Mpc and (0.48±0.10, 0.30±0.07) for R > 10h” 1 
Mpc, and hence are not significantly sensitive to the min¬ 
imum scale used. The \ 2 of the best-fitting constant Eg 
model to R > 10 h -1 Mpc is (5.5, 5.6) in the two redshift 
bins, for 6 degrees of freedom. The consistency of these 
measurements with the predictions of the perturbed GR 
model demonstrates that this model has successfully pre¬ 
dicted both the scale-independence of the signal, and the 
overall amplitude. 

Our measurement agrees with the fit obtained by Reyes 
et al. for z = 0.32, Eg = 0.39 ± 0.07. The errors obtained 
in the value of Eg are also similar; the Reyes et al. mea¬ 
surements are based on SDSS imaging data which covers 
much wider area, but provides a significantly lower density 
of measured shapes. 

We now test the sensitivity of our measurements to two 
important analysis choices. First, in Figure fTTl we show the 
results of repeating these fits of Eg to the range R > 10 h -1 
Mpc, in the two redshift bins, for different choices of the 
parameter Ro used for small-scale suppression when mea¬ 
suring T g m and T gg , in the range 1 < i?o < 3 ft -1 Mpc. 
These results show that our measurement of Eg is insensi¬ 
tive to the choice of Ro. Secondly, in Figure [T51 we compare 
the fiducial measurements to repeated fits applying differ¬ 
ent cuts in the range of photometric redshifts included in 
the analysis: 0.2 < zb < 2.0, 0.4 < zb < 2.0, 0 < zb < 1-6, 
0 < zb < 1-2 and 0.2 < zb < 1-3. We also show results 
just including RCSLenS or CFHTLenS regions, and for a 
“blue” sub-sample selected by a cut Tb > 2 (which will be 
subject to a lower level of intrinsic alignments, as discussed 
in Section [A4I) . The measured values of Eg never differ by 
a significant amount from the fiducial case. 

Finally, we assess the significance of any bias in 
the mean or standard deviation of Eg induced by non- 
Gaussianity in its posterior probability distribution. Since it 
is a ratio of noisy quantities, the error distribution of Eg is 
not necessarily Gaussian, even if the numerator and denomi¬ 
nator are Gaussian-distributed. We compared the mean and 
standard deviation obtained from direct error propagation 
with the distributions arising from Monte-Carlo sampling 
the relevant quantities in the numerator and denominator. 
We found that, averaging over bins with R > 10/i _1 Mpc, 
the mean of Eg was biased low by A Eg = (0.03, 0.01) in the 
two redshift ranges 2 = (0.32, 0.57), which is just under half 
the statistical error for the lower-redshift bin (the greater 
impact in the lower-redshift bin is driven by the greater frac¬ 
tional error in /3). The width of the l-cr and 2-<r confidence 
regions agreed within 5% for the two cases. We conclude that 
systematic errors due to skewness in the Eg distribution are 
sub-dominant to statistical errors. 

An additional issue would arise if the numerator and 
denominator of the Eg relation are additionally correlated, 
since then the mean (or median or mode) of the ratio would 
not be an unbiased estimator of the true value (e.g. Viola et 
al. 2014, Marian et al. 2015). We have estimated the correla¬ 
tion strength between the galaxy clustering and the galaxy- 
galaxy lensing correlation functions using Gaussian analytic 
covariances and found the correlation to be negligible. This 
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Figure 16. Covariance matrix of the Eq(R) measurements in the 
two redshift bins, displayed as a correlation matrix C%j / ■sjC%%Cjj. 
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Figure IT. The dependence of the fit of constant Eg to the 
range of scales R > 10 h ~ 1 Mpc on the value of Ro chosen for the 
small-scale suppression in the annular differential surface density 
statistics. 


is driven by two effects: (a) the samples studied in this work 
are dominated by shot noise rather than cosmic variance, 
which tends to de-correlate two-point statistics, and (b) w p 
and AE are related to the underlying matter power spec¬ 
trum via different weighting functions, such that the same 
separations in w p and AE are sensitive to different ranges 

of Pmm(k). 


7 CONCLUSIONS 

In this study we have performed a new test of gravitational 
physics based on a consistency check between the amplitude 
of peculiar velocities, measured via redshift-space distortion 
in galaxy samples from the WiggleZ and BOSS surveys, and 
the galaxy-galaxy lensing signal imprinted in the shapes 
of background galaxies in the CFHTLenS and RCSLenS 
datasets by density fluctuations traced by these lenses. Our 
results agree with the predictions of GR for a perturbed 
FRW metric, in a flat ACDM Universe with matter density 
consistent with observations of the CMB. 
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Figure 18. The dependence of the fit of constant Eq to the range 
of scales R > ! 0 h 1 Mpc on the range of photometric redshifts 
allowed in the analysis. 


In particular we produce new measurements of the 
“gravitational slip” statistic Eg = T gm /(/3T gg ), estimat¬ 
ing values 0.48±0.10 at 2 = 0.32 and 0.30 ±0.07 at 2 = 0.57 
when averaging over scales 10 < R < 50 /i _1 Mpc, compared 
to model predictions of Eg = 0.41 and 0.36, respectively. 
The results are consistent with those obtained when aver¬ 
aging over 3 < R < 50/U 1 Mpc. When carrying out this 
measurement we suppressed small-scale information from 
R < R 0 = 1.5 h _1 Mpc using the annular differential sur¬ 
face density statistic, although we find that our results are 
in fact insensitive to the choice of Ro. 

Assuming a ACDM cosmological model and fixing all 
the parameters apart from as, we determine <78(2 = 0.32) = 
0.75 ± 0.08 and < t 8 (2 = 0.57) = 0.54 ± 0.07 by fitting to 
the differential surface density AS and projected lens corre¬ 
lation function w p , after marginalizing over separate linear 
bias factors of each lens sample. These results are also con¬ 
sistent with the expected growth of structure in a perturbed 
ACDM Universe. 

In terms of methodology, we particularly note the sys¬ 
tematic bias that can occur in the measurement of AE (R) 
when using photo -2 information to select source-lens pairs, 
and how this can be reduced by using the full redshift prob¬ 
ability of each source in the estimator. We find that our re¬ 
sults are not sensitive to the range of photometric redshifts 
included in the analysis. 

We note that a measurement of Eg using lensing of the 
Cosmic Microwave Background by BOSS-CMASS galaxies 
was recently presented by Pullen et al. (2015), determin¬ 
ing E g {z = 0.57) = 0.243 ± 0.060(stat) ± 0.013(sys). This 
measurement is consistent with ours, with both results lying 
below the standard model prediction. 

Combinations of gravitational lensing and redshift- 
space distortion data offer a powerful opportunity to test 
large-scale gravitational physics in the search for an under¬ 
standing of the physical nature of dark energy. Datasets such 
as CFHTLenS, RCSLenS, the Kilo-Degree Survey (KiDS), 
Dark Energy Survey (DES) and the Hyper-Suprime Cam 
(HSC) lensing survey, promise that the constraining power 
of such tests will continue to improve rapidly and provide 
increasingly stringent tests of fundamental cosmology. 
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Table 2. Values of constant Eq fit to the measurements for each combination of source-lens datasets, using the range of scales R > 3 h 1 
Mpc and R > 10/i —1 Mpc. Best-fitting values of x 2 are shown for each case, and the measurements of the input RSD parameter 3. 


Survey 

0 

Eo(R>3h x Mpc) 

X 2 /dof 

E G (R >10 h J Mpc) 

X 2 /dof 

WGZLoZ-RCSLenS 

0.65 ±0.10 

0.36 ±0.15 

16.5/11 

0.52 ± 0.21 

13.3/6 

LOWZ-RCSLenS 

0.38 ±0.11 

0.40 ±0.14 

13.6/11 

0.54 ±0.15 

4.1/6 

LOWZ-CFHTLenS 

0.38 ±0.11 

0.41 ±0.16 

8.4/11 

0.36 ± 0.18 

4.9/6 

WGZHiZ-RCSLenS 

0.63 ±0.08 

0.38 ±0.15 

10.0/11 

0.27 ±0.17 

8.1/6 

CMASS-RCSLenS 

0.36 ±0.06 

0.36 ±0.08 

18.9/11 

0.36 ± 0.09 

13.3/6 

CMASS-CFHTLenS 

0.36 ±0.06 

0.24 ±0.09 

10.8/11 

0.23 ±0.11 

5.9/6 

Combined 0.15 < 2 : < 0.43 

- 

0.40 ± 0.09 

10.3/11 

0.48 ± 0.10 

5.5/6 

Combined 0.43 < 2 : < 0.70 

- 

0.31 ±0.06 

8.1/11 

0.30 ±0.07 

5.6/6 
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APPENDIX A: CALIBRATION CORRECTIONS 
AND SYSTEMATICS TESTS 

Al Systematics tests manipulating source shears 

We performed a series of systematics tests manipulating the 
source and lens catalogues and re-measuring the average 
tangential shear 7 1 ( 0 ): 

• We rotated the sources by 45° (i.e., ei, n ew = e 2 ,oid and 

^2,new = 6l ? old)* 

• We randomized the shears (i.e., randomly shifted each 
pair of values (ei,e 2 ) to a different source galaxy). 

• We replaced the lens catalogue by a random catalogue. 

The results of these tests are shown in Figure IaTI in all 
cases we should find 7 1{9) = 0 in the absence of systematic 
errors. The covariance matrices for each measurement are 
obtained by applying the same systematics-test operation to 
each mock catalogue; this covariance allows us to evaluate 
the \ 2 statistic of the -yt = 0 model in each case. These 
values are printed in each panel along with the number of 
degrees of freedom, dof = 20 , and indicate consistency with 
the model. 

A2 Shear bias calibration corrections 

The effect of the additive shear calibration bias is negligi¬ 
ble for galaxy-galaxy lensing. Figure IA2I displays the mul¬ 
tiplicative shear bias corrections K(9) and K(R), defined 
by Equations [33] and 1351 which are applied to the galaxy- 
galaxy lensing measurements by multiplying the estimated 
values by (1 +A) -1 . Corrections are shown for each combi¬ 
nation of source-lens datasets, combining the different sur¬ 
vey fields. These corrections are approximately independent 
of scale and have values K ~ —0.06. We note that the de¬ 
terminations of K(R) depend on the redsliift distribution 
of the source and lens samples, through the weighting fac¬ 
tor E^T 1 , which imprints an extra variation between datasets 
compared to K(9). 

A3 Photo-a systematic test: SIS fits in zb slices 

Systematic errors in the BPZ photometric redshift probabil¬ 
ity distributions would imprint errors in the determination 
of AE (R), which relies on the computation of E7 1 for each 
source-lens pair from an integral over redshift weighted by 
Pbpz(z), and in the cosmological modelling, which uses the 
source redshift probability distribution p s (z) derived from 
stacking the individual Pbpz{z) functions. We tested for 
such systematics by computing the amplitude of the tan¬ 
gential shear around a fixed set of foreground lenses in a 
series of eight photometric redshift bins defined using zb - 
(0 - 0.2,0.2 - 0.4, 0.4 - 0.6, 0.6 - 0.8, 0.8 - 1.0,1.0 - 1.2,1.2 - 
1 . 6 ) - in order to determine whether this amplitude scaled 
with redshift in the expected manner. A convenient method 
for quantifying the results is to assume a singular isother¬ 
mal sphere (SIS) model for the lenses, characterized by a 
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Figure Al. Systematics tests of the 7 1 ( 6 ) measurement applied to the galaxy-galaxy lensing cross-correlations for the different com¬ 
binations of source-lens datasets. The X-axis is the angular separation 9 in degrees, the y-axis plots 10 5 7 t• The left, middle, and right 
columns respectively show measurements following rotation of sources by 45°, randomization of shear values, and replacement of the 
lens catalogue by a random catalogue, x 2 statistics are quoted for the measurements with respect to a model of zero, indicated by the 
horizontal line. 


velocity dispersion a v , and to verify that the values of a v 
derived from the shear of each photo-z source slice are con¬ 
sistent. The shear profile for SIS lenses is given in terms of 
the Einstein radius 6e by 7 1(9) = 6e/26 , where 


9e — 4"7t 



and the geometrical factor is given by 


(Al) 
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jf-) = I dzpi(z 


dz p a (z') 


X(z') - \(z) 

X{z') 


• (A2) 


Figure [A3] displays the stacked RCSLenS BPZ redshift prob¬ 
ability distributions p s (z), weighted by the lensfit weights, 
in each of the photo-z slices. A signal is produced even for 
low values of zb, because sources are scattered to apparent 
low redshifts from true higher redshifts. 

Figure [ATI displays the fits for the a v parameter of the 
four different lens samples - WGZLoZ, WGZHiZ, LOWZ 
and CMASS - in each photo-z bin, using these p(z 3 ) func¬ 
tions determined from the outputs of BPZ. The fits are per¬ 
formed to 7 1 ( 9 ) measurements in the range 0.01° < 9 < 1°. 
The inferred values of <y„ are generally consistent as the 
source photo-z slice changes, with the possible exception of 
the lowest redshift slice 0 < zb < 0.2 for CFHTLenS, which 
contains a negligible number of galaxies. 


A4 Intrinsic alignments systematic test: SIS fits 
to Tb samples 

Another potential systematic effect that could impact our 
measurements is the intrinsic alignment of the background 
source sample with the foreground lens distribution. This 
effect is expected to preferentially diminish the average tan¬ 
gential shear of red source galaxies compared to blue galax¬ 
ies. The BPZ pipeline returns a source galaxy type in the 
range 0 < Tb < 6, and we used this parameter to divide 
the sources into sub-samples of red (Tb < 2) and blue 
(Tb > 2) galaxies. Following Section IA3I we then deter¬ 
mined the velocity dispersion of foreground lenses in the 
SIS model for the blue, red and combined source samples 
in the range 0.2 < zb < 1.6. We also performed mea¬ 
surements for a blue sample over a photo-z range signif¬ 
icantly higher than the maximum redshift of the spectro¬ 
scopic lenses (0.9 < zb < 1.6 for WGZHiZ and CMASS, and 
0.6 < zb < 1.6 for WGZLoZ and LOWZ. The source distri¬ 
butions Ps(z) for each case were determined by stacking the 
appropriate subset of BPZ redshift distributions. Results are 
displayed in Figure IX5l The inferred values of a v are again 
generally consistent amongst the samples, with the possible 
exception of the cross-correlations between CFHTLenS and 
BOSS red galaxies. However, as shown in Figure HB1 the fit¬ 
ted value of Eg is not significantly affected by cutting the 
sources to a blue sample. 


A5 Source-lens association correction B(R) 


Some sources may be clustered or associated with the lenses, 
but scattered to higher redshifts by photo-z errors. These 
sources will not be lensed, diluting the cross-correlation sig¬ 
nal. The strength of this effect may be determined by mea¬ 
suring the excess in the number counts of source galaxies 
in the vicinity of lens galaxies, compared to a random dis¬ 
tribution of lenses. The resulting multiplicative bias in the 
measurement of A'E(R) may be corrected by boosting the 
signal by 


B(R) = 


y y 

—^sources x ' 


data lenses j W i W j E c,^ ©(*. j) 


y y 

£—^sources i £—j 


sources i random lenses j 




(A3) 


where w r denotes the weights of the random lenses (nor¬ 
malized such that Y2j w j — w j)’ an< ^ where should 
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Figure A2. Multiplicative shear bias corrections to be applied 
to 7 /, (0) (upper panel) and A T,(R) (lower panel), for the different 
combinations of source-lens datasets. 
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Figure A3. Stacked RCSLenS BPZ redshift probability distri¬ 
butions, weighted by the Zerasfit weights, in a series of zg slices. 


be replaced by when moving from the estimator of 

Equation [30] to 1281 

Figure lA6l displays the boost factors B(R) for shapes 
in the different RCSLenS and CFHTLenS regions corre¬ 
lated with WiggleZ and BOSS lenses, averaging over 40 
random catalogues. The signal is generally consistent with 
B = 1, with deviations at the level of 3% or below for scales 
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R > 2h~ 1 Mpc. We do not apply this correction to our 
measurements. 

We note that an effect B(R) < 1, which is not expected 
through physical association, may be obtained if there is 
an anti-correlation between sources and lenses (for exam¬ 
ple, star-forming galaxies avoid dense environments on small 
scales), or through instrumental effects impacting the imag¬ 
ing data (such as contamination from light of foreground 
galaxies or sky subtraction). 


-200 0 200 400 600 

ct v [km/s] 



Figure A4. The SIS velocity dispersion fit to the average tan¬ 
gential shear around a series of foreground lens samples, for dif¬ 
ferent sets of sources in photometric redshift slices split by zb- 
The source redshift distributions for each slice, needed to model 
the resulting signal, are obtained from the stacked BPZ redshift 
probability distributions obtained for each source. 



Figure A5. The SIS velocity dispersion fit to the average tangen¬ 
tial shear around a series of foreground lens samples, for different 
sets of sources split by galaxy type according to T#. The source 
redshift distributions for each slice, needed to model the result¬ 
ing signal, are obtained from the stacked BPZ redshift probability 
distributions obtained for each source. 
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Figure A6. Boost factor B(R) from source-lens clustering for the different combinations of source and lens samples analyzed in this 
study. The measurements are the average over 40 random catalogues, with the errors indicating the standard deviation between the 
catalogues. 
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